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Abstract 

Topology on the lattice is reviewed. In quenched QCD topological 
susceptibility \ is fully understood. The Witten-Veneziano mechanism 
for the 7/ mass is confirmed. The topological susceptibility drops to 
zero at the deconhning phase transition. Preliminary results are also 
presented for \ an d x' i n fuh QCD, and for the spin content of the 
proton. The only problem there is the difficulty of the usual Hybrid 
Monte Carlo algorithm to bring topology to equilibrium. 



1 Introduction 

Topology plays a fundamental role in QCD. The key equation is the Ua(1) 
anomaly 

<V* = -2jv>g(x) (i) 

where = Yli=i ^Cf'l^i is the singlet axial current, Nf the number of light 
flavours and 

9 2 

Q( x ) = Q^F^Fp^^po- ( 2 ) 

is the topological charge density. Q(x) is related to the Chern current K^(x) 

d.K^x) = Q(x) (3) 

with 

K, = (dp^ - \gf ahc A b p Al} (4) 
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and as a consequence Q = f Q(x) is an integer on smooth classical configu- 
rations with finite action. Eqs. (0-0) yield 

d,(4(x) + 2N f K,(x))=Q (5) 

whence Ward identities can be derived. At the leading order in the 1/N C 
expansion, (N c — ► oo, with g 2 Nf fixed), Q(x) is zero, Ua{X) is a symmetry 
and rf is its Goldstone particle, = 0. The anomaly acts as a perturbation 
and shifts the position of the pole from zero to the actual rf mass. From 
Eq.(|5|) by inserting the quark mass terms, the relation follows [[[], 

where x is the topological susceptibility of the unperturbed vacuum (1/N C = 
0) defined as 

X = J d\x - y)d*dy(0\T(K,(x)K u (y))\0) quenched . (7) 

The subscript "quenched" indicates that the matrix element has to be com- 
puted on the ground state of the 1/N C = theory. In particular this implies 
that fermion loops, which are 0(g 2 Nf) are put to zero. Eq.(|5]) gives 

X = (180 MeV) 4 (8) 

which is expected to be valid within an order 0(l/N c ) of accuracy. Eq.@ is 
a peculiar equation relating physical quantities (masses, f n , Nf) to x, which 
exists in an artificial 1/N C = world. Its verification, however, is a check 
of the validity of the expansion, which is a fundamental issue. Lattice is 
an ideal tool to produce this artificial world, and in particular the absence 
of fermions in it simplifies the numerical work. Eq.(|7|) uniquely fixes the 
prescription for the singularity in the product of operators K ll (x)K u (y) as 
x — > y: 5-like singularities disappear after integration and this uniquely 
determines x- Eq.fl/]) is a specific prescription for the notation 

X = J d 4 x(0|T(Q(x)Q(0))|0> 

quenched • (9) 

When computing x by any regularization scheme, like lattice, an appropriate 
subtraction must be performed to satisfy the prescription Eq.(|7]). 
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The behaviour of \ at finite temperature, and more specifically, at decon- 
finement, is an important key to understand the structure of QCD vacuum || 

A regularized version of the operator Q(x), Ql(x), can be defined on 
the lattice. There is a large arbitrariness in this definition, by terms of 
higher order in the lattice-spacing which go to zero in the continuum limit. 
In general Ql{x) will not be an exact divergence, so that a multiplicative 
renormalization Z with respect to continuum can exist If a is the lattice- 
spacing, 

Q L (x) = a 4 ZQ(x)+0(a 6 ). (10) 
A lattice susceptibility xl is defined as 

Xl = ^(Ql(x)Q l (0)) = @P-. (11) 

x V 

In general the definition ([LI]) will not satisfy the prescription Eq. ([/]). It will 
be[l 

XL = a 4 Z 2 x + M + 0(a 6 ) (12) 

where M is a mixing with the continuum operators having dimension < 4 (/3 
is the beta function) 

M(J3) = A(D)iM-F'^)a l + P(I3) ■ 1. (13) 



X = ^ (14) 



From Eq.Qg 

V = — 

Xl is measured on the lattice numerically, a 4 is determined as usual by com- 
parison to a physical quantity (p mass, string tension); M and Z can be 
determined non-perturbatively by a procedure known as heating |J. The 
idea is that classical configurations with known topological charge, can be 
dressed by quantum fluctuations without modifying the topological content 
since topological charge is difficult to change by the usual local Monte Carlo 
algorithms. Z is determined by measuring the total charge Ql on a config- 
uration with an instanton, where Q — 1, Eq. fllOl) . M is determined from 
Eq. (|T2|) by measuring xl on the sector Q = 0, where as a consequence of 
Eq.©, X = 0. 
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Figure 1: x f° r SU(3). Diamonds, circles and squares correspond to the 0, 1 
and 2-smeared operators. 
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Figure 2: x f° r SU(2). Diamonds, circles and squares correspond to the 0, 1 
and 2-smeared operators. 



In Eq.fll4T) xl, M and Z strongly depend on the choice of Ql, on the 
action and on the coupling constant (3 = 2N c /g 2 . a depends on the action 
and on (3. x must be independent of all these parameters. Fig. 1 shows x 
for SU(3) 0, determined for several /3's and by use of different operators: 
as visible in the figure % = (175 ± 5 MeV) 4 . For SU(2) (Fig. 2) [| X is 
somewhat larger: x = (198 ± 6 MeV) 4 . 

The determination by using the so-called geometrical method, if accom- 
pained by the appropriate subtraction, agrees with the other choices of Ql- 

Fig. 3 shows the behaviour of x across deconfinement. The drop is 
stronger for 377(3) than for SU{2) @. 

2 Full QCD 

By use of the same procedure described in the previous section, x can b e 
determined in full QCD. The expectation from the Ward identities is that 

jv> _ 

X~ (Y^m^i) +0(m 2 q ). (15) 
i=i 
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Figure 3: \ across the deconfinement transition for SU (2) (circles) and SU (3) 
(squares). 

Our preliminary result, extracted from simulating at (5 — 5.35 where a = 
0.11(1) fm with 4 staggered fermions at am = 0.01 is 

X = (110 ±8 MeV) 4 (16) 

to be compared to the predicted value 

771 — 

X=^(^) m =o~(109MeV) 4 . (17) 

On the same sample of configurations we obtain for x' the preliminary value 

x' = 258 ± 100 MeV 2 or ^ = 19 ± 4 MeV (18) 

which is compatible with the value expected from sum rules || a/x 7 = 25 ± 
3 MeV. However both these determinations are preliminary because of the 
effect shown in Fig. 4 where we display the history of the topological charge 
Q along the Monte Carlo updating which produces the configurations []T0f| . 
In the updating algorithms used in quenched QCD (Metropolis, heat-bath) 
the topological charge has tipically 50 — 100 steps of authocorrelation time. It 
thermalizes slowly with respect to local quantum fluctuations (and this is the 
basic property which allows the heating method for the measurement of Z and 
M as explained in section 1), but a thermalized sample of configurations can 
be prepared in a reasonable CPU time. The algorithm used with dynamical 
fermions, the hybrid Monte Carlo, performs very badly in that respect, as is 
visible from Fig. 4. The configurations there correspond to about 700 CPU 
hours of an APE Quadrics with 25 GFlop which is a huge time. Our sample 
has thus a much smaller number of independent configurations than shown 
in Fig. 4, and therefore the errors in the results given in Eqs. (|I~6"D and (|I~8|) 
are underestimated. 
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Figure 4: History of Q in a Hybrid Monte Carlo run as a function of the 
molecular dynamic time r. 



The same uncertainty affects our determination, on the same sample, of 
the spin content of the proton. The matrix element of between proton 
states can be parametrized as 

(p*, s'\J% s) = utf, s>) (G^k 2 )^ 5 - G 2 (k 2 )^%) u{p, s) (19) 

with k = p — p'. The form factor G\ is related to the so-called spin content of 
the proton AE, Gi(0) = AE where AE = Au+Ad+As is the contribution of 
the different quarks species to the spin of the proton. The naive expectation 
would be AE ~ 0.7. The value determined from the moments of the spin 
dependent structure functions of inelastic scattering of leptons on nucleons 
is much lower: AE = 0.2(1). The lattice allows a determination of AE from 
first principles. One possible technique consists in the direct measurement 
of the matrix element (|T9|) . An alternative is to use the anomaly equation, 
which after taking the divergence of both sides of Eq.©, 

(i?,s'\Q\p,s) = ^AtfWifrfirfu&s) (20) 

A(k 2 ) = G 1 (k 2 ) + —G 2 (k 2 ). (21) 

m N 

As k — > Eq. (PT|) determines Gi(0), unless G 2 (fc 2 ) has a pole at k 2 = and 
this is the case in the quenched approximation but not in full QCD. Eq. ([2~0|) 
gives thus AE in terms of the matrix element (p 1 , s'\Q\p, s), which can be 
measured on the lattice. In principle the lattice operator Ql would mix with 
d^J^x) and but this mixing, as well as the small anomalous dimension 



of Q can be neglected JTT|. Our preliminary value is AE = 0.04(4). Here 
again the error could be larger and in any case the value is preliminary, due 
to the bad sampling of topology in our ensemble of configurations. 
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3 Conclusions 



Measurement of the topological susceptibility x on the lattice is fully under 
control. For quenched SU(3) the value is in good agreement with the pre- 
diction of I], |2|. x drops to zero at the deconfining transition. Preliminary 
determinations of x' m full QCD agree with sum rules. The spin content of 
the proton is at hand. The practical problem is the thermalization of topol- 
ogy on the lattice. Our huge sample of configurations is not thermalized with 
respect to it. This creates in principle a problem for the lattice determination 
of any quantity: a priori, indeed, it is not known how it could depend on 
the topological sector and therefore if the ensemble is biased with respect to 
topology, this could affect the result in an impredictable way. Solutions of 
this numerical problem are currently under study. 
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